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ABSTRACT 



Supersonic flow past oscillating cylindrical shells and 
oscillating flat-plate cascades is analyzed using linearized 
characteristics methods. Pressure distributions and gener- 
alized aerodynamic forces are computed for arbitrary radius 
to length ratios, axial and circumferential wave numbers, 
Mach numbers, and reduced frequencies. Excellent agreement 
is obtained with a previous solution of the complete un- 
steady linearized potential equation by Dowell and Widnall 
using Laplace transform techniques. Also, pressure dis- 
tributions are computed for flat-plate cascades with sub- 
sonic leading edge locus for arbitrary solidity, stagger 
angle, frequency, and interblade phase angle. For compar- 
ison, a two-blade solution is developed using the method 
of singularities. Good agreement is obtained for this 
special case. 
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I. 



INTRODUCTION 



Supersonic flutter of skin panels first occurred on the 
German V-2 missile and has become of increasing engineering 
importance in recent years with the development of super- 
sonic flight vehicles. Panel flutter, for example, was 
found to occur on the X-15 Rocket Aircraft and Saturn Space 
Launch Vehicle. An integral part of the flutter analysis of 
such panels is the determination of the oscillatory pressure 
distribution. A number of approaches have been formulated 
which resulted mostly in approximation theories. A detailed 
review of the state of the art can be found in Reference 6 . 
Only one "exact" approach was worked out by Dowell and Widnall 
using the complete linearized unsteady potential equation 
(Ref. 2) . However, this approach requires rather complicated 
and time-consuming complex integration techniques. 

In the present work a method of characteristics procedure 
is developed to determine supersonic flow past vibrating 
panels and cylindrical shells resulting in a rapid computation 
of the pressure distribution and generalized aerodynamic 
forces . 

A second flutter problem has also become of considerable 
interest very recently. Due to the continuing trend toward 
higher tip speeds in modern jet engine fans and compressors, 
the possibility of supersonic blade flutter has become of 
considerable concern. A recent review of this problem area 
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..■Oi;, 



« 



m 



has been given in "A VJorkshop on Aeroelasticity in Turbo- 
machines" (Ref. 13). 

In the present study the method of characteristics is 
applied to compute supersonic flow past a finite flat plate 
cascade with subsonic leading edge locus. In the past, cas- 
cade studies have proved to be a useful tool for the under- 
standing of the complicated three-dimensional, incompressible 
or subsonic flow through an actual compressor stage. Al- 
though it is doubtful that the much more complex mixed 
transonic-supersonic flov/ can be modeled by a two-dimensional 
cascade idealization, the analysis of an infinitely thin flat 
plate cascade may point the way toward a fuller understanding. 
The method of characteristics again proved to be a convenient 
tool for the analysis of supersonic flow past oscillating 
cascades . 
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II. 



PROBLEM FORMULATION 



Inviscid adiabatic flow of a gas along the positive x- 
axis is considered which is perturbed by either A) an air- 
foil cascade executing small pitch oscillations (Fig. 1) or 
B) small amplitude vibrations of the surface of a circular 
cylinder whose axis is aligned with the free-stream direction 
(Fig. 2) . 

The following basic equations hold: 

The equation of state for a perfect gas , 



p = pRT 

the conservation equation for mass, 

Ei + p (il + il + !T) = 0 



Dt 3x 9y 

the conservation of momentum 



9 2 



and the conservation of energy 



DU ^ 


1 


= 


Dt 


P 


9x 




1 




Dt 


P 


9y 


Dw 


1 


= 


Dt 


P 


92 


energy 




Ds 


0 





( 1 ) 

(2) 



Dt 

D 

where Dt is the substantial derivative and s is the entropy 
of the fluid particle. 

The assumption of small amplitude vibrations makes it 
possible to introduce the small perturbation concept. 
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Figure 2. Cylinder Geometry 
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Therefore, all flow quantities are assumed as small pertur- 
bations linearly superimposed on free-stream quantities, 

i.e. , 

u = + u (5) 

V = V (6) 

W = w (7) 

c = c„ + c (8) 

and the pressure and density perturbations are; 

Ap = P-Poo (9) 

Ap = p-p„ (10) 

The pressure perturbation can be expressed by the velocity 
of sound perturbation, i.e., 

_ 2 /-X 

P*“Poo Y — 1 (C""Coo) 



( 11 ) 



Specializing to planar flow, one has for the continuity 
equation; 



2 3c . 2 3c , 3u , 3v 

+ u- — c„ — + — =0 



( 12 ) 



Y”1 3t Y~1 3X 

and the Euler equations ; 



3x 



3y 



3U 


+ 

oo 


!^ + 


2 

z-' 


3c 


3t 


8x 


^ oo 

Y-1 


3x 


3v 


+ 




2 

n 


3c 


3t 


3X 


8 

) 

1 

>- 


3y 



(13a) 



0 (13b) 

Equations (12) and (13) are the governing equations for the 
case of supersonic flow past oscillating cascades. 

For the case of supersonic flow past vibrating cylindrical 
shells, the equivalent single equation for the perturbation 
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I 



potential, <I>(x,r,0,t) in cylindrical coordinates was taken 
as the more convenient basic governing equation, i.e.. 



(1-M2)$ +$ +1$ $ -2- $ = 0 

' ^ vv rr r r r'^ 00 c xt c*^ tt 



XX 



(14) 



For a derivation of this equation reference is made to one of 
the standard texts, e.g., Bisplinghof, Ashley, and Halfman, 
"Aeroelasticity" , (Ref. 1). 

For a complete formulation, proper initial and boundary 
conditions must be prescribed. The most important ones are 
A) the flow on the body surface must be tangent at any instant 
in tim.e and B) Sommerfeld's radiation condition must be satis- 
fied, i.e., waves must propagate away from sources of 
disturbance. 

After properly linearizing the general flow tangency 
condition, (see Ref. 1) one obtains for the flow tangency 
condition on a thin airfoil, performing a time dependent small 
amplitude motion about the y=0 plane. 



3h^ 3h 

V = — — 

3t 3X 



at y=0 (15) 



where h(x,t) is the equation of the airfoil chord line. 

For the vibrating cylindrical shell the flow tangency 
condition over the flexible portion becomes 



$ 

r 




r=R 3t 3x 



0<x<l 



$ 



r=R 



0 



x<0 



(16) 
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where h(x,9,t) in this case describes the radial shell 
deflection . 

With the assumption of a perturbation potential, the 
momentum equations, equations (3) , may be integrated to 



3<J>,u 3$.rdp 

+ °° +u^ + — -=c (t) 

3t 2 ~3x i p 



( 17 ) 



(18) 



assuming uniform parallel flow far upstream, 

f.. u^ 

C (t) =— 22. 

2 

Also, the density-pressure relationship can be expanded as 

2 , 



P (P)=P„+ (P-Pa>)+I (|^) = (P-Pco) ^+ - • • 



( 19 ) 



therefore , 



■P dp ^1 


•P 


-| _P~Poo 


•? p (p) p 


p 


_ P 



oo oo 



dp: 



.P-Pa 



( 20 ) 



Hence, the perturbation pressure can be expressed in terms 
of the perturbation potential, i.e.. 



p-p = -p 



3$j 3$ 

— +U 



(21) 



3t 3x 

In pressure coefficient form this becomes 



p r 



TP 

^ 00 00 



u5 t X 



( 22 ) 



For further analysis it is convenient to use non-dimen- 
sional quantities. Introducing for the oscillating cascade 
case the blade chord as the unit of length one has the 
following non-dimensional coordinates: 
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X - y r 

X = — , V = — , t = 



t u. 



( 23 ) 



c c c 

Furthermore, introducing the assumption of harmonic time- 
dependence, one can define a reduced frequency 

(24) 



k = — 



where w is the circular frequency of blade oscillation. 

Following previous v;ork by Teipel (Ref. 3) for single 
blade oscillations the non-dimensional complex amplitude 
functions 

T /.4+- 11 

(25) 



TT/- 

u (x,y)e 



u 

u,. 



icx)t 

V(S,y)e 



^ iwt 
C (x,y)e 



V 



/ M ^-1 

2 1 c 

Y~1 Co 



(26) 



(27) 



are introduced. Substituting these equations into the continu- 
ity and Euler equations (12) and (13) gives 



|£ + / m ^-1 |X + -t- ikM^C = 0 

ax ay ax 

« + |£ + iku = 0 
ax ax 



av 



1^ + ikv = 0 



(28) 

(29) 

(30) 



v/M^ri 9y 

where the overbar on x and y is being dropped for convenience 
in the subsequent analysis. 

Introducing similar non-dimensional variables with dis- 
tances referred to flexible cylinder length, velocities to 
free-stream speed and time to cylinder length divided by 
free-stream speed, one obtains the equation for the non-dimen- 
sional perturbation potential, i.e.. 
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J 



0 



(31) 



(1_M2)$ + $ + 1 $ _ 2m2$ - 

XX rr 99 xt tt 

subject to the boundary condition 

<I>„I 



r=R 



ah 3h 
at ax 



0<x<l 
= 0 x<0 

If the radial shell deflection is assumed as 

ilc t 

h(x,e,t) = Z (x) cos n0*e 
where k is now referred to cylinder length, i.e. 

wL 



(32) 



k = 



u. 



(33) 



(34) 



Then the velocity potential must also have the form, 

ilc "h 

$(x,r,6,t) = (})(x,r) cos n0*e (35) 

Substitution of equation (35) into equation (31) then gives 

(1-m 2)4, + ^ + ^ - 54 ({, - 2ikM2^, + = 0 (36) 



subject to 






rr 



r r 



r=R 



ax 



+ ikZ 

= 0 



X 



0<x<l 

x<0 



(37) 



The pressure coefficient in terms of the non-dimensional 
velocity potential is 

C = -2 [ik$ + $ ] (38) 

p ^ X 



and the pressure coefficient amplitude is 



Cp = -2[ik({> + 4> ] (38a) 

These are the basic equations used in the subsequent cylind- 
rical shell vibration analysis. 
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III. 



THE METHOD OF CHARACTERISTICS 



A. CYLINDRICAL SHELL 

In Section II, equation (36) was developed as the partial 
differential equation for <j>(x,r) that requires solution 
subject to the boundary condition, equation (37) . The follow- 
ing developments are based on a previous paper by Platzer, 

Brix, and Webster, (Ref. 7). 

If M>1, equation (30) is hyperbolic and has real character- 
istics, which satisfy the ordinary differential equation 

(M2-l)dr2-dx2 = 0 (39) 

An arbitrary function F(x,r) has the following derivatives 
along the two characteristics: 



^ = F + F dr 

ds ^x*ds ^r ds M 

where ds is the differential 

istic ds2 = M^dr^ and 



• F ± ^ F (40) 

X M r 

arc-length along the character- 



dx _ / M^-1 
ds ~ M 



(41) 



and 



dr _ + 1 
ds “ M 



(42) 



Letting s ^ be the arc- length of the left-running Mach line 
and 52 the arc -length of the right-running Mach line, one 
obtains for the first derivatives. 



Fi = 



/ m -^-1 „ 

M X 




(43) 
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and 



F 2 = -^TT^ F - i F 

M X M r 



and for the second derivatives 



( 44 ) 




Figure 3. Cylindrical Shell Characteristics 

Network 



The finite difference form of this equation for the points 
P and S (see Figure 3) on the right-running Mach line is 

't’llS) - (P) ^ ^ ^ 



As 



2rM /M^-1 



(<j>l+<{>2) + (^^ 4> (47) 



similarly, for the points Q and S on the left-running Mach 
line one obtains 

ikM 



^2 (S) ~ 4'? (Q) _ <{>i“<i>? /, ^ , /,. 

Ai 2m /m^ (4>i+<f>2) + (k 



2__H1 



.2m2 






(48) 
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Now solving for the values of <{>1 and ((>2 at S from the known 
values at P and Q (Figure 3) the entire flow field may be 
computed point by point. For this purpose the values for <j)i 



and (j )2 in equation (47) are averaged between P and S, i.e., 

(49) 

. (P) + 4>2 (S) 

4,2 - 2 (50) 



. _ <^1 (P) + (S) 

^ 



and 



4>(S) = (}.(P) + ^[<}. 2 (P) + <J. 2 (S)]As (51) 

is obtained by integration along the right-running Mach line. 

Similarly, <1) j and equation (48) are averaged be- 

tween Q and S 



^ 1 (Q) + <1> 1 (S) 

(f)l - — I 1 



(52) 



^ ? (Q) + 4' 2 (S) 

<P2 “ 9 



(53) 



and 

4>(S) = (J)(Q) + i [(|)i(Q) + 4.1 (S)] As ' (54) 

2 



This results in a system of two equations for 4.1 (S) and 4.2 (S) , 
i.e . , 

A4>i (S) + B4>2 (S) = C (55) 

D 4.1 (S) + E4>2 (S) = F (56) 
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where 



^ _ T As , ikMAs 

^-mrcsT 2m^ 



B =- 



As ikMAs _ As2 2 J N 

^msr " 2/fpri „2,2(s)’ 

^ n . As ikMAs, 

- <l> 1 ( ) [ + 4Mr (P) “ 2 /m^ 

. ^ iT^\ r As ikMAs . As^,, o 

^2(P)[- 4MFW- - 



(57) 



(58) 






-)] 



D = 



As 



ikMAs As^ ^^,2 



n‘ 



2/M^-l M^r^CS) 

As ikllAs 



E = - 1 - 



®FTsT 2 /h^ 

A s . ikMA s - 



M^r^ (S) 
(59) 

•)] (60) 
(61) 



F = <^2 (Q) [-1 + ~T^ + _^_^] 

™(Q) 2/M^ 

. . r„ AS , ikMAs A 

wrm* ”4 



■(k2 - 



n‘ 



+ 4>(Q) [- 4^(k2 - 









(S)M^ 



) - |2.{k2 - 

^ **2^2 



M^r^ (S) 
2 



-)] 



-)] 



M^r2 (Q) 

(62) 

The above procedure will suffice for the computation of 
(|» and its derivatives at a general flow field point, but needs 
to be modified to obtain the values on the initial Mach line 
and along the cylinder surface. 

On the initial Mach line (Figure 3) it is known that 
({,=({, j=0. Hence, equation (56) becomes 

E<^2(S) = F (63) 
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where E is known from equation (61) and F simplifies to 



F 



<}> 2 (Q) 



[-1 + 



As 

4Mr (Q) 



+ 



ikMAs ^ 

2v^M^-l 



(64) 



On the cylinder surface the normal velocity is prescribed, 



i.e. , 

<l>r = I = If + 

whereas equation (55) must hold along the right-running Mach 
line, The values for <J>i and ({>2 on the surface are then found 
from 



4>i (S) 



* I 'H + 

A + B 



( 66 ) 



and 



C - 



<f>2 (S) = 






A + B 



(67) 



The pressure distribution over the vibrating cylindrical 
shell is obtained from the linearized pressure amplitude, 
equation (38a) , and the generalized aerodynamic force 0^^^ ^ is 
found from 

Q = 2 / C (x) Y (x)dx (68) 

rnr 2 0 Pjj^ t 

where 



Z(x) = Z A (x) 
mm 
m 



(69) 



and Cp^(x) is the aerodynamic pressure amplitude due to the 

deflection 'i' (x) . 

m 

These equations were programmed in Fortran IV. The 
program description appears in Section IV-A. 
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B. CASCADE 



In Section ii equations (28) - (30) were developed as the 
governing system of partial differential equations for the 
analysis of the cascade problem. A solution of these equations 
can be obtained using the method of characteristics. The 
following development of this method is largely based on the 
previous work by Teipel (Ref. 3) and Chalkley (Ref. 9), v;here 
further details can be found. 

Transforming the system into a new coordinate system, 

C(x,y) and n(x,y), and considering now a system of three simul- 
taneous equations in unknov/ns U^, , and C^ , the equations 

may be written 

c + cot a K Cj, = - U„ny “ “ ^n^y 

X 5 y ? X ? n ^ 

- M^c n - ikM^C (70) 

n X 

CU, + CC, = -Un - Cn - ikU (71) 

X C X ? n X n X 

C + tan a^C, =-Vn - tan a C n - ikV (72) 

X c y C n X n y 

If U^, , and C^ are made indeterminate across lines of 

n = constant, the determinant of the matrix of coefficients 
of the above equations must equal zero, or 

- - C [C ^ - m2 ^2] = 0 (73) 

X X y X 

The solutions to this equation give the three character- 
istic directions in the physical (x,y) plane: 




^ = ± tan a (75) 
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In order to find the compatibility relations of the first 
derivatives of U, V, and C along the characteristics, equations 
(70) - (72) are written along linos of 5 = constant# n = 
constant, and along streamlines, that is, for 



C(x,y) = constant, 
n(x,y) = constant. 



tant , 


0 

9x 


(76) 


= 


1 


(77) 


9x 


/m^-1 


= 


- 1 


(78) 


9x 


/m^-1 



If an arbitrary function of x and y in the equations is de- 

9 f 

noted by f then, t — holding str constant becomes 

o X 



(i£) 

9x str 



9X 



(79) 



9f 



Likewise, — holding C constant becomes, 

o X 






3f 



9x ^ 9x /[ 4 2 _ 2 ^ 



(80) 



9f 



and, holding n constant becomes 

o X 



( 1 ^.) = ~ = 



9f 



n 9x 

Combining these equations and rewriting. 



9f _ . 9f > 

9x 9x^str 



(81) 



— = ^ [ ( 
9x T ' 



il.) + (M) ] 

9x'C 9x%^ 



9y 



n—). - (~) ] 

2 - 9x C 9x n 



,9f, 



(82) 

(83) 

(84) 



Substituting these relations into the equations of motion 
gives 
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1 f, 3 U. ^ ,au, , , 1 ^ 

2 ‘Wc ^ <n'n’ * 2 '<35?’ 5 ■ <35?’ n’ 



+ 7 m2 [ (|£) 

2 '3x 


5 * 




0 


(85) 


1 [ {H) + (i!i) ] 

2 ^^9x^c ^9x n 


+ 


1 [ (^) + (1£) ] 
2 ^'9x^c ^9x n 


+ iku = 0 


' (86) 


t(^) + (^) ] + 

^^9x^5 ^9x n 


1 

2 


[ (i£) - (i^) ] + 

^^9x^5 ^9x n 


ikV = 0 


(87) 



Subtracting equation (86) from equation (85) gives an equation 
that is first added to equation (87) and then subtracted from 
equation (87) . The result is two compatibility relations that 
must be satisfied along lines of constant C and lines of 
constant r\r i.e., 

( 88 ) 

(89) 



^) + 
9x^ C 


(^) 

'9x'C 


+ ik[V + — — 
m2-1 


(m2c- 


-U)] = 


^) - 
9x' n 


(9£) 

^9x^ n 


+ ik[V - — — 
m2-1 


(m2c- 


-u)] = 



The third compatibility relation is obtained by considering 
equation (86) along a streamline 



(|^) +. + (t^) +. + ikU = 0 

9x str 9x str 



(90) 



If irrotational flow is assumed, the equation of irrota- 
tionality 



9u _ 9v _ Q 
9y 9x 



(91) 



may be substituted for the second equation of motion in the 
system. ' With the assumption of simple harmonic motion, it 
may be v/ritten in terms of Teipel's amplitude functions. 



9y 9x 



(92) 



The non-dimensional system of equations can then be written 
as , 
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7 . 



+/ + m 2 + ikM2c 

3x ay ax 



+ ikU = 0 

ax ax 



|ii = 0 

ay ax 



(93) 

(94) 

(95) 



By transforming equations (93) - (95) to the arbitrary co- 
ordinate system, C and n / the system can be considered one of 
simultaneous equations in U^, , and • By making any of 

these derivatives indeterminate, the results of setting the 
determinant of the matrix of their coefficients equal to zero 
are equations (74) and (75), the characteristic directions. 

By writing the equations of motion, as before, along 
lines of C = constant and n = constant and performing the 
same algebraic manipulations, the two compatibility relations 
can be determined: 



, au > 

ax^ C 



IV) 
ax' 5 


+ ik 

m2-1 


(U-C) = 0 


(96) 


^)' 
ax' n 


, m2 

m2-1 


(U-C) = 0 


(97) 



, au . 

ax^ n 



The third compatibility relation is, as before, 

Ostr = 0 <90) 



The equations of motion are nov>7 reduced to a system 
containing only derivatives with respect to x. It is, thus, 
possible to solve these equations using finite differences in 
the following manner: separating the real from the imaginary 

parts, the equations can be written. 
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0 



(98) 



9x ' 


+ 




str 


aUj 




ac 


Tx”^ 


+ 


^9x 


str 



Rn 



R 



9U, 



9V. 



3U. 



9V, 



2 . 

-). + k (U - C ) = 0 

5 m2-1 ^ ^ 

T) m2 

^) - k (U^ - c_) = 0 

^ m2_i I I 

I 

i) + k (U^ - C^) = 0 

n i42_i R R 



( ■*■) 
'ax ' 5 


— 


( 

'ax 


( — 2.) 
'ax ' n 




av. 


+ 


<9X 


au 




av. 


‘sx >n 


+ 


( 

'ax 



(99) 

( 100 ) 

( 101 ) 

( 102 ) 

\ 

(103) 



To V7rite these equations in finite difference form the 
computational molecule shown in Figure 4 v;ill be used. With 
the velocity values all known at P 2 . 1 ' ^12' ^21' 

values at ^22 calculated. The distance Ax is arbitrary. 




Figure 4. General Flow Field Molecule 
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If F is an arbitrary flow quantity, its partial deriv- 



atives can be written as 



F - F 

.8F. • _ 22 11 

'Dx^str Ax 


(104) 


;9Fv _ ^22 " ^12 

C 1 * 

Ax 


(105) 


,9F. ^22 " ^21 

I-4X 


(106) 



The values of the flow quantities in the equations are 
averages , such that 



(I')jtr 2 *'22* 


(107) 


(F)^ = i (Fi2 + F 22 ) 


(108) 


‘'■'n = ? ‘''21 ^ '' 22 > 


(109) 


With these relations substituted into equations 


(98) - (103), 



the system becomes after some algebraic manipulation. 



^22R ^22R “ ‘^1^221 ” ^12R 


(110) 


^221 ^221 ^I^22R ^121 


(111) ' 


^22R " ^22R ®I^^22I " ^22l^ " ^34R 


(112) 


^221 " ’^221 ®I^^22R " ^22R^ ^ ^341 


(113) 


^22R ^22R ®I^*^22I “ ^221^ ” ^56R 


(114) 


^221 ^221 ®I^^22R " ^22R^ ^ ^561 


(115) 



where 



^12R ”i 1R ^IIR ^l^lll 


(116) 


^121 " ^111 ^111 " •^I^llR 


(117) 
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(118) 



34R 


^12R 




^12R 


+ 




<“121 


" ^121 


341 


^121 


- 


^21 


- 




<^2R 


^12R 


56R "" 


^21R 


+ 


^21R 


+ 




*'^211 


^211 


561 


^211 


+ 


^211 


- 




*'^21R 


^21R 



and 



Aj = j kAx 
1 



= 5- k 



M‘ 



Ax 



(119) 

( 120 ) 
( 121 ) 

(122) 

(123) 



m2-1 

As shown in Teipel (Ref. 3), solving these equations 
gives , 

^22R "" ^F^^34r'^^ 56R^ ”^1^121^ ■'■^®I^2‘^^34l'''^56I^‘‘‘^I^12R^ 



(1-AjBj)2 + (2Bj)2 



(124) 



U 



221 



(1-AiBi) (^34 i'''^56I^'' ®i^12R^"^®I^2 ^^34r''’^56R^ "®I^12I^ 



(1-AjBj)2 + (2Bj)2 

^22R = 7 (^6R - ^34r) 



^221 2 ^^561 ” ^341^ 



^22R ^12R ” ^22R ^1^221 



^221 ^121 " ^221 " '^ I ^ 22 R 



(125) 

(126) 

(127) 

(128) 
(129) 



These are the finite difference equations for a point in the 
general flow field. 

At a point on the upper surface of a cascade blade, the 
normal velocity "^22 prescribed by the movement of the blade. 
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This velocity was given in the problem formulation. Here the 
computational molecule is shown in Figure 5 . 




X 



Figure 5. Upper Surface Molecule 
The applicable equations are, 

^22R *^22R " ■^1^221 ^ ^12R 

^221 ^221 ^l’^22R ^121 

^22R “ ^1^221 ^1*^221 ^ ^56R " ^34R 

^221 ®I^22R " ®I^22R "" ^561 " ^341 

where _ and K^, are as before, but, 

±Z t)0 

= \7 

34R 22R 



^341 ^221 



(130) 

(131) 

(132) 

(133) 



(134) 

(135) 



as given by the boundary conditions. Solving for the flow 
quantities gives 



^22R ^^56r"^34r“®I^^ 12I^'’'^®I ^^56l"^34l'*'®I^12R] 



( 1 - AjB ^)^ + ( 2 B ^)2 



(136) 



U 221 - (l-AA)[K5gj-K3^^+B3.K^2Rl-2Bi[K3g^-K3^j^-BjK^2ll 



(1-AiBi)2 + (2Bj)2 



*^22R ^12R “ ^22R ^1^221 



(137) 

(138) 
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(139) 



^221 ^121 " ^221 “ ^ I ^ 22 R 



and ^ 221 ^ ^221 known from the boundary conditions. 

These are the finite difference equations for the flow 
velocities on the upper surface of a cascade blade. 

Likewise, at the lower surface of a cascade blade, the 
normal velocity, '^22 prescribed by the movement of the 
blade, and was given in the problem formulation. The compu- 
tational molecule, different from that at the upper surface 
of the blade, is shown in Figure 6. 




Figure 6. Lower Surface Molecule 
Here, the applicable equations are. 



^22R ^22R “ ■^1^221 ^12R 



^221 ''' ^221 ^ I ^ 22 R ^121 



^22R ■ ?I^22I ®I^22I ^56R ^34R 

^221 ®I^22R " ®I^22R ^ ^561 ^341 



where K ^2 are as originally stated, but. 
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^56R = ^22R 
^561 "" ^221 



(140) 

(141) 

(142) 

(143) 

(144) 

(145) 
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as given by the boundary conditions. Thus, the finite 
difference equations for the flow velocities at the lower 
surface of a cascade blade are, 



U 



^l-Ml)t^56R+^34R-Vl2I^-^2B3.[K3g3.+K3^3.-fB^K^2R] 



22R 



U 



221 





(1-AjBj)2 + (2Bj)2 


(146) 


(1-Aj.Bj) tKjgj 


■'■^34l‘''®I^12R^ f^56R’‘'^34R“ 


®I^12I^ 




(1-AjB^)2 + (2Bj)2 


(147) 


c 

^22R 


^12R “ ^22R ^1^221 


(148) 


c 

^221 


= ^121 - ^221 - V 22 R 


(149) 



where V 22 j^ snd '^ 221 . given from the boundary condition. 

To obtain the conditions on the initial left-running 
Mach line, assume that and ^21 shown in Figure 4 are 

just in the free-stream and then let Ax shrink to zero. 

Since K ^2 ^^55 = 0 ^^id B^ go to zero as Ax goes to 

zero, the initial finite difference equations for C=constant 
are , 



thus , 



“22 <^22 = « 



“22 - '^22 * ''12 = 0 



“22 '^22 ® 



U = - V = - c 

22 ^22 ^22 



(150) 

(151) 

(152) 

(153) 



Equation (96) can be written, 

<l^>, = - 



(154) 
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Integrating gives. 






U = U| »e”^^M2-l 
' x =0 



•X 



thus , 



= -V„__(0)cos (k 



M' 



22R 



22R 



M 2-1 



x) -V 221 ( 0 ) sin (k 






m 2-1 



(155) 



X) (156) 



U 



221 



= -V__t( 0 )cos (k 



M' 



221 



M^-l 



x) ( 0 ) sin (k 






22R 



m 2-1 



^22R ^22R 



^221 ' ^221 



C = - u 
^22R 22R 



C = - U 
^221 ^221 



x) (157) 

(158) 

(159) 

(160) 
(161) 



For the initial right-running Mach line, and P ^2 

Figure 4 are assumed to be just in the free-stream and Ax 
is brought to zero. Since K ^2 “ ^34 = 

(162) 
(163) 



“22 <^22 = 0 



thus , 



“22 - ''22 = “ 



“22 - “21 * '^22 = ° 



^^22 ^22 ~ *'22 



Equation (97) can then be written. 






Integrating gives , 



U = V(0)e 



.•1, 

-ik 

m 2-1 



(164) 



(165) 



(166) 



(167) 



34 



Thus , 



“22R = ''22K<0>‘^°= "" 



u 



221 



V22i(0)cos (k 



MZ 


x) +V 221 (0) 


(k 


m2 


x) 


(168) 


M2-1 


m2-1 


m2 




(k 


m2 


X) 


(169) 


m2-1 


m2-1 


22R 


^22R 








(170) 


22R 


“^22R 








(171) 


221 


^221 








(172) 


221 


"’^221 








(173) 



At the first position encountered on a blade other than 
the first blade there are two grid points for which the flow 
field quantities must be computed; one on the upper surface 
and one on the lower surface of the blade. 

As in Figure 7 the blade is assumed to protrude a very 
small distance past P 22 ' necessitating the computation 

of the upper and lower flow velocities at this point. In this 
case, the equations developed earlier for an upper surface 
grid point are used for F 22 U those developed for a lower 

surface grid point are used at 
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At a position just aft of a blade there are also two 
grid points to be accounted for; one above and one below the 
slip plane in the wake. As in Figure 8 the point 
assumed to be a very small distance aft of the blade's trail- 
ing edge. 




Figure 8. Wake Molecule 



At this position the equations for an upper blade grid 
point apply to ^2211 those for a lower grid point apply to 

^22L' 



22RU 


+ 0 — A n = 

22RU I 22IU 


^12RU 


(174) 


22IU 


+ r + A II = 

22IU I 22RU 


^12IU 


(175) 


22RU 


®I^^22IU " ^22IU^ 


^56R “ ^22RU 


(176) 


22IU 


®I^^22RU " ^22RU^ 


- K - V 

561 22IU 


(177) 



and 



22RL 


4- 


*^22RL " ■^I^22IL 


^12RL 




(178) 


22IL 


+ 


^22IL ^I^22RL " 


^12IL 




(179) 


22RL 


+ 


^*“22IL " ^22IL^ 


^22RL 


^34R 


(180) 


22IL 




®I^^22RL “ *“22RL^ 


^22IL 


^341 


(181) 
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where 



12RU 


^IIRU 


‘^IIRU 




(182) 


12IU 


^llIU 


^llIU 


■^1^1 IRU 


(183) 


12RL 


= ^IIRL 


^IIRL 


VllIL 


(184) 


12IL 


= ^IIIL 


■*' ^IIIL 


- ^I^llRL 


(185) 


34R 


^12R - 


^12R ^ 


1^^121 “ ^121^ 


(186) 


341 


^121 - 


^121 - ® 


I^^12R “ *^12R^ 


(187) 


56R ^ 


^21R 


^21R ® 


1^^211 “ ^211^ 


(188) 


56l 


^211 


^211 - ® 


I^^21R - ^21R^ 


(189) 



Equations (174) - (181) form a system of eight equations v/ith 
twelve unknowns. However, since it is known that the wake 
cannot support a pressure difference and that the normal 
velocities on either side of the wake must be equivalent, i.e.. 



^22RU 


= V 

22RL 


^ ^22R 


(190) 


^22IU 


^22IL 


^ ^221 


(191) 


^22RU 


^22RL 


^ ^22R 


(192) 


^22IU 


^22IL 


" ^221 


(193) 



we can eliminate four of the unknowns and achieve a determinate 
system of equations. Equations (174) - (181) become 



^22RU 


+ 


^22R 




A II = 

I 22IU 


^12RU 


(194) 


^22IU 


+ 


^221 


+ 


A II — 

22RU 


^12IU 


(195) 


^22RU 


+ 


^22R 


+ 


®I<=22I - 


^22IU^ “ ^56R 


(196) 


^22IU 


+ 


^221 


+ 


^I^^22RU 


^22R^ ^561 


(197) 
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(198) 



^22RL ^22R " •^I^22IL ^12RL 



^22IL 


*^221 


H- 


^I^22RL ~ ^12IL 


(199) 


^22RL 


^^22R 


+ 


^1^^221 ” ^22IL^ ^ ^34R 


(200) 


^22IL 


^221 


+ 


^I^^22RL " ^22R^ ^ ^341 


(201) 



The pressure distributions along the surfaces of the blades 
can be determined from equation (11) . Dividing both sides of 
this equation by and using Teipel's amplitude function for 
c, the equation can be written, 

P (K,y) = YM2c(x,y) (202) 

where 

P(x,y)e^*^^ = ^ (203) 

P« 

Thus the non-dimensional pressure distributions along the 
upper and lower surfaces of the blades in a cascade can be 
computed in a point by point calculation of U, V, and C along 
lines of C=constant and n=constant. 

These equations were programmed in Fortran IV. The pro- 
gram description appears in Section IV-B. 



C. THE LINEARIZED METHOD OF CHARACTERISTICS OF SAUER AND HEINZ 
Sauer and Heinz (Ref. 11) developed a linearized method 
of characteristics for supersonic flow past bodies of revor 
lution at zero angle of attack. This method also is derived 
from the linearized potential equation for steady supersonic 
flow, i.e., 



(m2-1) ^ 



XX 




(204) 
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Sauer and Heinz showed that the compatibility relations can 
be written as 



i^)= r ||- (205) 

and 

iilZL „ /fprr r (206) 

3n 9n 



where 




(207) 



V 



ii 

ar 



(208) 



and 

C = X “ /m^- 1 r (209) 

n = X + /M^-i r (210) 

are the characteristic coordinates as explained previously. 

Writing this equation in finite difference form one 
obtains for u and v at point 4 (Figure 9) using the known 
values in points 2 and 3 , - . 



^4 = 



(rv) (ry) (r^U 2 +r^U 2 



/m^- 1 (r_+r ) 
C n 



( 211 ) 



r_ (rv) -+r (rv) ..-/m^-I r_r (u_-u-) 
_ _J 3 n 2 g n 2 3 

^ r^ + r 

C n 



( 212 ) 



where r^ and r^ are the average r coordinates between 3 and 
4 and between 2 and 4 respectively. Computing the flow pro- 
perties at the body using the exact flow tangency condition 

% 

leads to 
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(213) 



^5 = 



/h 2-1 r^u^+(rv)^-rj (§|) ^ 



,dh. 



/M^-l (gj) 5 



''5 “ <il>5<^+“5> 



(214) 



where r = h(x) describes the body contour. 

However, in order to apply the Sauer-Heinz procedure to 
the case of an infinitely long cylinder v;ith a wavy surface 
between 0<x<l a different starting procedure is needed. It 
is easily seen from Figure 9 that 



r^ = “O (. 



2 /M-^-l tan 6 



+ 1 ) 



(215) 



0 



realizing that 



where 



and 



one obtains 



r = 
n 



^1 + ^2 



^1 = ^0 ^1 



^2 ^0 ^0 



r = R + ^ (3 + 



/M^-1 tan 6 



) 



(216) 



(217) 



(218) 



(219) 



0 



Since there is no disturbance upstream "'of the leading edge • 
Mach line, one obtains from equation (213) 



U. 



- (R+hp) tan 6 q 



/M^-1 r + (R+h ) tan 6 



( 220 ) 



where the point 2 is the first body point. 
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Finally, the pressure coefficient is obtained from 

C - 2u - (M^-1) (221) 

P 



The program listing for these equations is given in 
Appendix I. 
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IV. 



COMPUTATIONAL PROCEDURE 



A. CYLINDRICAL SHELL COMPUTER PROGRAM 



This program calculates the pressure distribution over a 
harmonically vibrating cylindrical shell. It is valid only 
for a cylinder in a supersonic flow field vibrating with 
small amplitude vibrations. 

The program calculates the flow field properties adjacent 
to the cylinder using the method of characteristics finite 
difference equations developed previously in Section III. 

The pressure distribution is then integrated using the trape- 
zoidal rule according to equation (68) to obtain the general- 
ized aerodynamic force. 

The computational diagram shown in Figure 3 is used to 
calculate the flow field quantities at S in the general flow 
field and on the cylinder surface. The distance As shown in 
the figure is determined by 



As = 



L cos a 



( 222 ) 



where v is the grid fineness ratio. The grid fineness ratio 
is an arbitrary input variable into the program and is equal 
to one less than the total number of grid points on the 
cylinder surface. 

Input parameters are entered into the program on two data 
cards. The first card contains the date, e.g., 11 APR 1973, 
in the first twelve columns. The second card contains the 
following input data: 
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SPACE 



DATA 



1-20 


Free-stream Mach number, M 


21-30 


Radius to length ratio, R 


31-40 


Reduced frequency, k 


41-45 


Axial mode number , m 


46-50 


r, as defined in equation (68) 
for the generalized aero- 
dynamic force. 


51-55 


Circumferential mode number, n 


56-60 


Grid fineness ratio, v 



where the first three quantities are floating point numbers 
and the rest are integers. The definitions of these para- 
meters are given again in the alphabetical listing of all 
program variables in Appendix A. 

The program has five subroutines. They are as follows; 



COMPXR 


Computes the x and r position 
of the given grid point. 


MACHLN 


Computes the flow quantities 
along the initial left-running 
Mach line at the given grid 
point. 


GENFPT 


Computes the flow quantities at 
a general flow field grid point. 


PANPNT 


Computes the flow quantities on 
the cylinder surface at the given 
grid point. 


INTEG 


Computes the generalized aero- 
dynamic force. 



Flow quantities along the initial left-running Mach line 
are computed using the equations derived in the method of 
characteristics for initial conditions. Those quantities on 
the cylinder surface are computed using the equations derived 
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for the boundary condition. The values of the velocity 
potential and its derivatives are derived as a result of the 
analysis of the flow tangency condition and have been given 
previously. The procedure starts at the initial Mach line 
and continues down successive right-running Mach lines. The 
pressure coefficient is computed when the cylinder surface 
is met and stored in an array for later use in the integration 
of the generalized aerodynamic force. 

A flow diagram of this program is given in Appendix B. 

B. OSCILLATING CASCADE PROGRAM 

This program calculates the pressure distribution over 
the upper and lower surfaces of the blades of a subsonic 
leading edge locus cascade which is oscillating in pitch at 
arbitrary frequency and interblade phase angle. It is valid 
only for a cascade in a supersonic flow field vibrating with 
small amplitude vibrations. The program calculates the flow 
field quantities within the area governed by the aft-running 
Mach lines emanating from the leading edge of the first blade 
and the forward-running Mach lines emanating from the trailing 
edge of the last blade of the cascade (see Figure 10) . These 
flow field properties are calculated using the method of 
characteristics finite difference equations developed previously 
in Section III-B. 

The computational molecules shown in Figures 4, 5, 6, 7, 
and 8 are used to calculate the flow field properties in the 
general flow field, upper surfaces, lower surfaces, leading 
edge of a staggered blade, and the wake aft of a blade. 
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respectively. The distance Ax shown in the figures is 
determined by 

Ax = ^ cot a (223) 

where v is the grid fineness ratio. This parameter is an 
arbitrary input variable into the program and is equal to one 
less than the total number of grid points on any Mach line 




Input parameters are entered into the program by means of 
the NAMELIST option of Fortran. In using this option, input 
parameters are combined in a titled list and then referred 
to in the program by this title (NAMl) . The real advantage 
of using NAMELIST is that each data card need only contain 
the value of one input variable while the format of the data 
card is such that the variable's name is punched on the card 
as well as its value. In this way, when a different value 
for an input variable is desired, rather than changing the 
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entire data deck, only one easily identifiable data card need 
be changed. 

The format for a NAMELIST data deck is as follows: The 
first column in each card of the data deck is left blank. 

The first card starts in the second column v?ith the symbol & 
(ampersand) followed immediately by the title of the NAMELIST 
(here, NAMl) . The next cards list the variables and their 
input values, one per card, in any order. The format here, 
starting in the second column and follov/ed by a comma, is: 
variable name = value. The last card in the data deck and the 
signal to the computer that the NAMELIST has ended again 
starts in the second colum,n with the symbol & immediately 
followed by END. 

The value of an input variable may be written in any 
format as long as integer variables are written without deci- 
mal points and real variables are written with a decimal point. 
Blanks are taken as zeroes, but a comma must appear somewhere 
between a desired input value and the next data card. The 
date of the run is entered on a data card in the first twelve 
spaces. This data card is inserted prior to the NAMELIST 
dataset. 

To illustrate the input procedure, the following example 
is proposed: 



V = 100 


d = 0.8 


M = /2 


O 

ID 

in 

II 

oa 


k = 0.75 


6 = 45* 



5 blades 
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The data 


cards needed are as 


follows ; 


CARD 


FORMAT 


DESCRIPTION 


1 


22 JUNE 1973 


Date 


2 


&NAM1 


NAMELIST title 


3 


NGRDFN=100, 


Grid fineness ratio, 
must be an even integer, 
subject to the constraint 
of equation (227) . 


4 


FSTRMN=1. 414214 , 


Free-stream Mach number 


5 


PvEDFRQ=0.80, 


Reduced frequency 


6 


XSUBO=0.50, 


Elastic axis position 


7 


. TNWDST=0.80, 


Vertical distance between 
blades . 


8 


STGANG=55.0, 


Complimentary stagger angle, 
tan B > cot a . 


9 


FAZE=45.0, 


Interblade phase angle 


10 


NUMBLD=5 , 


Number of blades 


11 


&END 


End of NAMELIST 



The definitions of these parameters are given again in the 

alphabetical listing of all program variables in Appendix C. 

The program has nine subroutines. They are as follows: 

INPUT Reads in all input data and sets up the finite 

difference grid. 

INTIAL Initializes all the flow quantities at the leading 

edge of the first blade and initializes most of 
the logic variables. 

MACHLN Computes the flow quantities along the initial left 

and right-running Mach lines at the given grid 
point. 

TOP Computes the flov? quantities on the upper surface 

of a blade at the given grid point. 
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V 



GENFPT 


Computes the flow quantities at the given general 
flow field point. 


BOTTOM 


Computes the flow quantities on the lower surface 
of a blade at the given grid point. 


NEWBLD 


Computes the flow quantities at the leading edge 
of all staggered blades. 


WAKE 


Computes the flow quantities in the wake aft of 
each blade at the given grid point. This sub- 
routine uses IBM subroutine SIMQ for the solution 
of its simultaneous equations. 


COMPXY 


Computes the x and y position of the given grid 
point. 



In INPUT the grid is determined using the grid fineness 
ratio, blade spacing, number of blades, and the free-stream 
Mach number. The value of Ax as shown in equation (223) is 



obtained . 


. The stagger angle is then adjusted such that the 


distance 


the upper blades are staggered back is the nearest 


multiple 


number of Ax increments. The "compatible stagger 



angle" is then printed out. The blades may be staggered such 
that the first staggered blade is at least one Ax increment 
aft of the initial left-running Mach line, but no further 



forward , 


or, 

d • (tan 3 - cot a)>Ax (224) 


Flow 


quantities along the initial left and right-running 



Mach lines are computed using the equations derived in the 
method of characteristics, Section B, for initial conditions 
Those flow quantities along the blades' surfaces are computed 
using the equations derived for the appropriate boundary 
condition. The equations for the normal velocities are de- 
rived as a result of the analysis of the flow tangency con- 
dition and have been given previously. The flow quantities 
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in the wake aft of each blade are computed using the set of 
simultaneous equations derived as a result of the assumptions 
made about the flow quantities on either side of the v/ake and 
are given in the Method of C haracteristics , Section III-B. 

The procedure starts at the leading edge of the first 
blade and computes the flow quantities at grid points along 
the initial right-running Mach line. It then jumps back to 
the initial left-running Mach line and computes the flow 
quantities along the next right-running Mach line, using the 
applicable subroutines at the appropriate grid points . Only 
flow quantities on the Mach line presently being computed 
and the Mach line which preceeded it are kept in storage at 
any one time. Flow quantities at successive grid points are 
computed in like manner until the trailing edge of the last 
blade of the cascade is met. 

A flow diagram of this program is given in Appendix D. 

C. SAUER-HEINZ COMPUTER PROGRAM 

This program computes the pressure distribution over a 
cylindrical shell or body of revolution. It is written so 
that the shape of the body or cylindrical surface may be 
easily changed. It is valid only for non-oscillating bodies 
or shells in a supersonic flow field at zero angle of attack. 
The following discussion applies to the cylindrical shell. 

The program calculates the flow field adjacent the 
cylinder using the method of characteristics finite difference 
equations developed previously in Section III-C. 
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The computational diagram shown in Figure 9 is used to 
compute the flow field properties at 4 in the general flow 
field and at 5 on the cylinder surface. The distance As 
shown in the figure is determined by, 

cos a 



As = 



(225) 



where v is the grid fineness ratio. The grid fineness ratio 
is an arbitrary input variable into the program and is equal 
to one less than the number of grid points on the initial 
Mach line. The length As^, n=l,2,3... is determined by the 
intersection of the right-running Mach lines with the surface. 

Input parameters are entered into the program on one data 
card. This card contains the following input data: 

SPACES DATA 



1-10 

11-20 

21-30 

31-40 

41-50 



Free-stream Mach num±)er , M 
Maximum amplitude of the axial mode, h^ 
Cylinder radius to length ratio, R. 
Axial mode number , m 
Grid fineness ratio, v 



where the first three quantities are floating point numbers 
and the last two are integers. The definitions of these 
quantities are given again in the alphabetical listing of 
all program variables in Appendix E. 

The program has three subroutines. They are as follows: 



COMPXR 



GENFPT 



Computes the x and r position of the 
grid point and determines the position 
of the intersection of the current 
right-running Mach line with the body. 

Computes the flow quantities at a general 
flow field grid point. 
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SURFAC Computes the flow quantities at the 

given surface point. 

Flow quantities on the initial left-running Mach line 
are assumed to be zero. Those quantities on the surface are 
computed using the appropriate equations developed in 
Section III-C, for surface grid points. 

The procedure starts at the initial Mach line and con- 
tinues down successive right-running Mach lines. The pressure 
coefficient is computed when the surface is met. 

A flow diagram of this program is given in Appendix F. 
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V. 



RESULTS AND COMPARISONS 



A. CYLINDRICAL SHELL 

Pressure distributions and generalized aerodynamic forces 
were computed for sinusoidal deflection modes, i.e., 

= sin(miTx) m=l,2,... (226) 

Figures 11 through 16 show results for the generalized aero- 
dynamic forces and Q-^2 ^ function of circumferential 

mode number for tv;o radius to length ratios, r-0.5 and r=1.0, 
at a Mach number of /2 and reduced frequencies of zero and 
one. Also shown- is the comparison with two previous solutions 
by Dowell and Widnall (Ref. 2) and by Matsuzaki and Kobayashi 
(Ref. 5). The Dowell and Widnall solution is also based on 
the complete linearized unsteady potential equation, but uses 
Laplace transform techniques for its solution. As can be seen, 
agreement between this solution and the method of character- 
istics solution is excellent. For the sake of interest, a 
recent approximate solution by Matsuzaki and Kobayashi, valid 
only for large values of the circumferential mode number is 
also shown. 

Figures 17 through 24 show the generalized aerodynamic 
forces srid Q 2 ^ 2 r ^ function of radius to length ratio 

again for a Mach number of ^2, zero reduced frequency, and 
circumferential mode numbers n=0,l,2, and 3. The main purpose 
of these figures is to show that the infinite cylinder theory 
is a poor approximation for small radius to length ratios, but 

approaches the exact theory for large values of R . 
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Figure 11. Generalized Aerodynamic Force vs. Circumferential Mode Number 
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Figure 12. Generalized Aerodynamic Force vs. Circumferential Mode Number 
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Figure 13. Generalized Aerodynamic Force vs. Circumferential Mode Number 
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Figure 14. Generalized Aerodynamic Force vs. Circumferential Mode Number 
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Figure 15. Generalized Aerodynamic Force vs. Circumferential Mode Number 
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Figure 16. Generalized Aerodynamic Force vs. Circumferential Mode Number 




Figure 17. Generalized Aerodynamic Force vs. 

Radius to Length Ratio 



59 




Figure 187 Generalized Aerodynamic Force vs. 

Radius to Length Ratio 
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Figure 19. Generalized Aerodynamic Force vs. 

Radius to Length Ratio 
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Figure 20. Generalized Aerodynamic Force vs. 

Radius to Length Ratio 
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Figure 21. Generalized Aerodynamic Force vs. 

Radius to Length Ratio 
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Figure 22. Generalized Aerodynamic Force vs. 

Radius to Length Ratio 
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Figure 23, Generalized Aerodynamic Force vs. 

Radius to Length Ratio 
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Figure 24. Generalized Aerodynamic Force vs. 

Radius to Length Ratio 
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Finally, Figure 25 shows the change in pressure distri- 
bution as a function of radius to length ratio for the case 
M=/2, m=l, n=0, k=0. These distributions were computed using 
the characteristics method outlined in Section III-A and the 
Sauer-Heinz theory as described in Section III-C. Complete 
agreement was found between the two methods. 

B . CASCADE 

Figures 26 through 29 show the instantaneous pressure 
distributions on the first and second blades of a finite 
cascade for four different values of reduced frequencies. 

The pressure difference between lower and upper blade surface 
is plotted for the instant when the airfoil is at its mean 
horizontal position and pitching nose up. Results for these 
cases were given very recently by L. E. Snyder, of Pratt and 
Whitney Aircraft (Ref. 8), for the following parameter 
combinations : 

M = 1.5 

Solidity = 1.33 (d = 0.317) 

Stagger angle = 25° (3 = 65°) 

Xq = 0.5 

S = 40° (Blade 2 leads blade 1) 
k = 0.04,0.332,0.832,1.25 

It can be seen that appreciable differences exist in a 
number of cases, in particular in the wave reflection regions. 
These differences are not fully understood at this time. 
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Figure 25. Pressure Distribution Over Wavy-Walled Cylinder 
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Figure 26. Cascade Pressure Coefficient Difference 

vs. Fraction Chord 
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Figure 27. Cascade Pressure Coefficient Difference 

vs. Fraction Chord 
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Figure 28. Cascade Pressure Coefficient Difference 

vs. Fraction Chord 
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Figure 29. Cascade Pressure Coefficient Difference 

vs. Fraction Chord 
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Snyder uses a finite difference procedure and further work is 
required to evaluate the validity and accuracy of this pro- 
cedure is compared to the method of characteristics procedure 
developed in the present work. However, a check could be made 
for the first blade solutions for v;hich the method of singular- 
ities of Garrick and Rubinow (Ref. 12) is available. This 
solution was programmed according to the equations given in 
Appendix H and the results are shown in Figures 26 through 
29. The method of characteristics solution is in complete 
agreement vrith the Garrick-Rubinow prediction whereas consider- 
able deviations are seen to exist for the finite difference 
predictions . 

A further check using the method of singularities was 
developed for two blades oscillating with zero phase difference. 
This approach is also described in Appendix H and a comparison 
with the method of characteristics is given in Table 1 below: 



TABLE 1 



k 




•0 


.1 






0 


.5 




M 


/2 


1 


.8 


/2 


1 


.8 


X 


M.O.C. 


M.O.S. 


M.O.C. 


M.O.S. 


M.O.C. M 


.O.S. 


M.O.C. 


M.O.S. 


0.1 


.979 


.979 


.667 


.660 


.693 


.694 


.640 


.486 


0.2 


.977 


.977 


.666 


.661 


.686 


.687 


.638 


.511 


0.3 


.974 


.974 


.666 


.662 


.680 


.681 


.635 


.536 


0.4 


.972 


.972 


.666 


.663 


.676 


.677 


.631 


.562 


0.5 


.969 


.969 


.666 


.664 


.674 


.675 


.627 


.588 


0.6 


.966 


.966 


.665 


.666 


,675 


.676 


.623 


.614 


0.7 


.963 


.963 


.665 


.667 


.678 


.679 


.618 


.641 


0.8 


.959 


.959 


.665 


.668 


.682 


.683 


.614 


.668 


0.9 


.956 


.956 


.664 


.669 


.689 


.690 


.607 


.694 
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Further method of characteristics results for a finite 
subsonic leading edge cascade are shown in Figures 30 through 
32 for the follov;ing parameter combinations: 

M=/2 Xq=0.5 

d=1.0 k=0.5 

6=59.8° 6=120° 

As mentioned previously, the computer program is written 
sufficiently general so that the pressure distribution can 
be computed and printed out for arbitrary specified Mach 
number, reduced frequency, stagger angle, blade spacing, and 
number of blades subject to the following constraint: 

V tan a [ (N-1) (tan 6-cot a) + i]<800 (227) 

Equation (227) is the expression which governs the maximum 
grid fineness that may be used with a given set of input 
parameters. This constraint is necessary because the program 
is dimensioned such that the maximum number of grid points 
allowed on any one right-running Mach line is 400. 
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Figure 30. Pressure Coefficient Difference vs. Fraction Chord 
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Figure 31. Pressure Coefficient Difference vs. Fraction Chord 
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Figure 32. Pressure Coefficient Difference vs. Fraction Chord 



VI. 



CONCLUSIONS AND RECOMMENDATIONS 



A solution of the complete linearized unsteady potential 
equation has been developed for supersonic flov; past vibrating 
cylindrical shells using the method of characteristics. 
Pressure distributions and generalized aerodynamic forces 
were computed for shells of arbitrary radius-to-length ratios / 
axial and circumferential wave numbers, Mach number, and re- 
duced frequencies. A solution for this same problem was pre- 
viously obtained by Dowell and VJidnall using Laplace transform 
techniques. Excellent agreement was obtained in all cases. 

It is believed that the present method can be applied more 
easily than methods requiring complex variable techniques. 

It should be noted that the procedure can easily be modified 
to include arbitrary axial deflection modes. 

A further solution of the complete linearizjed unsteady 
potential equation using the method of characteristics was 
developed for supersonic flow past vibrating flat plate 
cascades with subsonic leading edge locus. Pressure distri- 
butions were computed for finite cascades of arbitrary 
stagger and interblade phase angles performing small ampli- 
tude pitch oscillations. Preliminary comparisons with a 
recent finite difference solution by Verdon and Snyder showed 
reasonable agreement. However, further work is necessary to 
assess the validity and accuracy of the present method. This 
should be done by further comparisons with the Verdon and 
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Snyder procedure and Fourier and Laplace transform methods 
currently under development at Stevens Institute of Technology, 
General Motors Corporation (Detroit Diesel, Allison Division), 
and General Electric. Also, flutter calculations should be 
carried out using the method of characteristics predictions 
of the oscillatory aerodynamic forces which would allow com- 
parisons with the Pratt and Whitney cascade flutter tests. 
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APPENDIX A 



CYLINDRICAL SHELL COMPUTER PROGRAM VARIABLES 
A. LOGIC VARIABLES 

The following variables appear in the logic statements 
of the program: 



IHAVEP 


The number of the computation along any 
right-running Mach line. IHAVEP is used as 
the first subscript for a grid point. 


ISWTCH 


Switch variable used to designate the 
right-running Mach line being calculated, 
alternately set equal to "1" or "2" and 
used as a second subscript. 


JSWTCH 


Switch variable used to designate the 
previously computed right-running Mach 
line, alv;ays the opposite of ISWTCH and 
used as a second subscript. 


KOUNT 


The value of IHAVEP at the cylinder surface. 


ML 


The number of grid points along the initial 
Mach line and along the cylinder surface. 



B. QUANTITY VARIABLES 

The following variables in the program take on the values 

defined below. All are floating point quantities unless 

otherwise stated. The dimensioned variables are listed first. 

AMAGCP(400) The imaginary part of the pressure coefficient 

CP (400) The complex pressure coefficient. 

PHI (400,2) The derivative of the velocity potential 

along the right-running Mach line. 

PHK(400,2) The derivative of the velocity potential 

along the left-running Mach line. 

QPHI (400,2) The velocity potential. 
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REALCP (400,2) 


The real part of the pressure coefficient. 


X,R(400 ,2) 


The X and r coordinates of a grid point. 


A,B,C,D,E,F 


See equations (57) - (62) 


BETA 


The square root of the quantity Mach number 
squared minus one. 


DELTAS 


The step size along the characteristics. 


DSTSTR 


The horizontal distance between grid points. 


FINGRD 


The grid fineness ratio. An even integer 
less than 400. 


HDSTRL 


One-half DSTSTR. 


HEIGHT 


The increment HDSTRL tan a . 


II 


The square root of minus one. 


IKM 

IKW 


The factor, ik — — - 

m2-i 

The quantity, ik sin mirx 


IPHIl, IPHI2 


The velocity potential as described in 
equations (51) and (54) . 


K 


The reduced frequency. 


MACH 


The free-stream Mach number 


MANGL 


The Mach angle, a. 


MFREQ 


The axial mode number. An integer. 


NFREQ 


The circumferential mode number. An integer 


NR 


r, as defined in equation (68). An integer. 


PARW 


The partial derivative of the deflection 
taken with respect to x. 


PI 


3.141593 


PX 


The partial derivative of the velocity 
potential taken with respect to x. 


RADIUS 


The radius to length ratio. 


W 


The cylinder deflection. 
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APPENDIX B 



CYLINDRICAL SHELL FLOW DIAGRAM 
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APPENDIX C 



OSCILLATING CASCADE COMPUTER PROGRAM VARIABLES 
A. LOGICAL VARIABLES 

The following variables appeal* in the logic statements of 
the program; 



ICO 


Set equal to "1" when a blade or V7ake com- 
putation is made. Set equal to "0" after 
the next general flow field computation 
occurs . 


IHAVEP 


The number of the computation along any 
right-running Mach line starting with zero 
at the initial left-running Mach line. The 
value IHAVEP+1 is the first subscript for 
a flov; quantity. 


ISOTCH 


Switch variable used to designate the 
right-running Mach line being calculated, 
alternately set equal to "1" or "2" and 
used as a second subscript. 


JLINE 


The Mach line counter incremented each time 
a right-running Mach line is completed. 
Reset equal to the value of IHAVEP at the 
leading edge of a new blade. 


JSWTCH 


Switch variable used to designate the 
previously computed right-running Mach 
line, always the opposite of ISWTCH and 
used as a second subscript. 


LCOUNT 


The value of IHAVEP at a surface or wake 
grid point. 


LIMIT 


The maximum value of IHAVEP. 


LIMW 


The value of LCOUNT at the trailing edge 
of a blade. 


MAXI 


The value of JLINE on a right-running Mach 
line v/hich intersects a blade's leading 
edge . 


NEWDST 


The number of Ax increments in a blade chord 
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NSTPTS 



The number of Ax increments between the 
initial Mach line and the leading edge 
of the first blade. 

NUM Set equal to one less than the number of 

the blade encountered on a Mach line. 

NUMBLD The number of blades in the cascade. 

NUMOLD Preserves the maximum value of NUM. 



NUMPI NUMOLD plus one. 

OLDL Preserves the value of LCOUNT at the highest 

blade or wake yet encountered. 



B. QUANTITY VARIABLES 

The following variables in the program take on the values 

defined below. All are floating point quantities unless 

otherwise stated. The dim.ensioned variables are listed first. 

A (8, 8) Matrix of coefficients in v;ake computation. 

B(8) Right-hand side of equations (194) - (201) 

U22R,U22I,V22R,V22I,C22R,C22I (400,2) . The velocities , real 

and imaginary of general flow field, upper 
blade, and upper v^ake points. Equivalent 
to the respective Teipel amplitude functions, 
^22R through ^ 221 ' 

U33R,U33I,V33R,V33I,C33R,C33I (400,2) . Same as above, but 

used only at lower blade and lower wake 
points . 

The X and y coordinates based on chord length. 

Aj- as defined in equation (122) . 

M^ 

The argument k ^x 

M^-l 

Bj as defined in equation (123) . 

The interblade phase angle in radians. 

Incremental distance between grid points 
along Mach line . 

The increment Ax. 



X,Y (400,2) 

AI 

ARC 

BI 

DELTA 

DELTAS 

DSTSTR 
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FAZE 



FMANGL 
FNGDPT 
FSTRI'IN 
Gl/ • • « /G9 

HDSTRL 

K12R, . . . ,K56I 
KS 

NGRDFN 

REDFRQ 

STAG 

STGR 

STGANG 

S 

TNVTOST 

TRNGLH 

U 

V 

VIPANL 

VRPANL 

W 



The interblade phase angle in degrees, 

Mach angle , a . 

Floating point NGRDFN. 

The free-stream Mach number. 

Intermediate factors used in velocity 
computations in GENFPT, TOP, BOTTOM, 
and NEWBLD. 



The increment Ax/2. 

^12R' ’ * ’ '^561 defined in equations (116) 
( 121 ) . 



Dummy variable used in SIMQ. An integer. 

The grid fineness ratio. An even integer, 
less than 400. 

The reduced frequency. 

The distance the second blade is staggered 
back from the y-axis. 

The distance from the initial Mach line to 
the second blade's leading edge. 

The stagger angle in degrees. 

The factor /m^-1 . 

The vertical distance betv/een adjacent blade 
The increment -^Ax cos a . 



The 


factor 


,, m2 
cos (k 


x) . 






m2-1 




The 


factor 


• /I m2 

sxn (k 


x) . 






•H 

1 

CM 





The imaginary part of the normal flov/ 
velocity at the first blade's leading edge. 

The real part of the normal flow velocity 
at the first blade's leading edge. 
m2 

The factor 

/ M ^-1 
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XSUBO 



The elastic axis position based on chord 
length. 



XLNGTH 

XN 

XNEW 



The length of the initial Mach line between 
the planes of adjacent blades. 

Floating point NUM. 

The x-coordinate on a blade based on the 
blade's leading edge. 
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APPENDIX D 



CASCADE FLOW DIAGRAM 
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APPENDIX E 



SAUER-HEINZ COMPUTER PROGRAM VARIABLES 



A. LOGIC VARIABLES 



The following variables appear in the logic statements 
of the program: 



I HAVE P 



ISWTCH 



JSWTCH 



KOUNT 



The number of the computation along any 
right-running Mach line starting with zero 
at the initial left-running Mach line. The 
value IHAVEP+1 is the first subscript for a 
flow quantity. 

Switch variable used to designate the 
right-running Mach line being calculated, 
alternately set equal to "1" or "2" and 
used as a second subscript. 

Switch variable used to designate the 
previously computed right-running Mach 
line , always the opposite of ISWTCH and 
used as a second subscript. 

The value of IHAVEP at the surface. 



B. QUANTITY VARIABLES 

The following variables in the program take on the values 
defined below. All are floating point quantities unless 
otherwise stated. The dimensioned variables are listed first 



CP(200) 

DELX,DELR(200) 

U(200,2) 

V(200,2) 

X,R(200,2) 

A,B 



The pressure coefficient. 

The incremental distances between grid 
points on a right-running characteristic. 

The streamwise velocity perturbation. 

The product vr. 

The X and r coordinates of the grid point. 
Used as factors in the equations for u and v 
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BETA 


The square root of the quantity M^-1. 


DELTAS 


The incremental distance betv;een adjacent 
grid points along the initial Mach line. 


DELTAX 


The change in x from one point to the 
next on the initial Mach line. 


DELTAY 


The change in r from one point to the next 
on the initial Mach line. 


FINGRD 


Floating point NGFN. 


HO 


The maximum amplitude of the deflection. 


MACH 


Free-stream Mach number. 


MANGL 


The Mach angle, a. 


MFREQ 


The axial mode number, an integer. 


NGFN 


The grid fineness ratio, an even integer, 
less than 200. 


NU14PTS 


The number of grid points on the initial 
Mach line. 


PI 


3.141593 


RETA 


r 

n 


RO 


The radius to length ratio. 


RZETA 




S 


The length of the initial Mach line. 


SLP 


Slope of the surface at a surface grid point 


XM 


Floating point MFREQ. 


XMAX 


The value of x at a surface grid point. 


YO 


The intercept of a right-running Mach line 
with the r-axis. 
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APPENDIX F 



SAUER-HEINZ FLOW DIAGRAM 
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APPENDIX G 



SUPERSONIC FLOW PAST AN INFINITELY LONG 
WAVY WALL CYLINDER 



Leonard and Hedgepeth (Ref. 4) gave a solution for super- 
sonic flow past an infinitely long wavy wall cylinder. De- 
tailed computations using this solvition were carried out by 
Anderson (Ref. 14). For the cylinder wall deflection 

h = hgSin m-JTX cos n 9 (G-1) 

where m is the number of half-waves per unit axial length 
and n indicates the number of circumferential waves, one 
obtains for the pressure distribution. 

Cp = 2h,. ^ A, sin(m7rx + (G-2) 

^ 0 M 1 1 



where 



, iY, M H 2 [mirR (m2-1) 2) 

Aj^e 1 = n 






1 



H 2- 2] 
n 



(G-3) 

For the actual computation equation (G-3) was rewritten as 

2M 



„ i’i'i 

Aj^e 1 = 



J ( Y ) “ i Y ( Y ) 
n ' n ' 



‘ Vi ' Vi Vi <''> > 



where 

Y = mTiR /M“^-l 

For the case n = 0 the equation was rewritten as 



(G-4) 

(G-5) 



i'i'T 

Aj^e 1 



-Jj^(y) + i Y^(y) 



M 



/M^-1 



(G-6) 
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and 



Ai = 



(G-7) 



T. = tan-1, 

Re(A^e"'*'l) 



(G-8) 



The program listing appears in Appendix I 
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APPENDIX H 



TWO BLADE SOLUTION USING THE 
METHOD OF SINGULARITIES 



As shown by Garrick and Rubinow (Ref. 12), the velocity 
potential of a single blade oscillating in supersonic flow is 






/ \ “1 ,x-y “ik^ 

(i(x,y)= — ■ ■ / v(5)e 

/ M ^-1 0 






kM 

M2-1 



V(X-C) ^-(M^-l)i^] dc 



(H-1) 




Figure H-1. Method of Singularities Configuration 



The solution for supersonic flow past two oscillating blades 
in Zone II, (Figure H-1), can be obtained by superposition of 
single blade solutions. For this purpose the disturbed flov; 
field in Zone II is written as the linear superposition of the 
single blade solutions for blades 1 and 2. However, in order 
to satisfy the flow tangency condition on the upper surface 
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} , 





of blade 2, a third single blade solution v/ith unknown normal 
velocity must be added. Hence, one can write 



<t>T "f*! *f'2 ‘^’iP 



(H-2) 



where 




r kM 

0 



V (X“C) ‘^-{M^-l)y‘^ ] d^ 



(H-3) 




= -1 ^Xi-y/M^'-l 



v(5)e (^1"5 ) (M^-1 ) yj^ ^ ] dc 



(H-4) 




(H-5) 



where denotes the unknown normal velocity. Equation (H-3) 
represents the disturbance due to the first blade, equation 
(H-4) the disturbance due to the second blade alone, and 
equation (H-5) the interference disturbance. Note that the 
blade 2 and interference potentials are expressed in terms of 
an x^, y^ coordinate system with origin at the leading edge 
of blade 2. 

Satisfying the boundary condition on the upper surface of 
blade 2 gives an equation for the unknown velocity, v^, i.e.. 



(x,y)= 






1 





/ (x-C) ^ - (M^-l)y"^] dC 



(H-6) 



at y = d . 
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Now setting 



V (O = ik Z (C) + Z' {a 



(H-7) 



and assuming the blades to pitch about their leading edges, 
i. e . , 

ZU) = K (H-8) 

and then performing the differentiation, one obtains 

. kM^ 



V (x,y)= — — ^ ^ ^ [ikZ (c)+ZM5) ] e 



kMy 



/ (x-^ ) (M^-1) y 



• J, /(x-C)^-(M^-l)y^]dc 

rruTT i M 1 



. kM^ 



- [ikZ (x-y/M^-1) +Z ' (x-y /M'^-l) ] • e ^M^-1^ 

at y = d (H-9) 

Since, in general 

= ik(}) + (J) (H-10) 

P X 

the total pressure coefficient can be written as 



S = Si S2 - s 



(H-11) 



IP 



where Cpj^ is obtained from the single blade solution for 
blade 1 at y = d. 






k 2 



•{/' 



/(x-o ^-{ M ^-l)d^]d?-i<|^/^"^‘^^^"^[-Z' {K)+^Z"iZ)] 



^0 
. kM^ 



(x-5) 






+ r-e M'^-l 
k 



(H-12) 
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( 

'I 



i 



( 






Cp 2 from the single blade solution for blade 2 at = 



1,2 .. . . 
" (/^l [-Z(C)+iz" (C) ]e 



o' 






kM^ 



(H-13) 



and 



Cp = ik $__ + ^ 



IP 



IP ^IP 



(H-14) 



where 

<j) 



-1 ,x 



kM^ 



/M"^-l 0 



: /^1 V^ (Y)e“^M^-l J^ [_|^(x^-y) ]dy 



(H-15) 



and 



-1 I ,x 
/M^-1 I 0 



f \ t • kM^ (x., -C ) 

/y vy y) [-ij^ZTT e ^ 



kM 

-r r kM - , -X- „ , . -l - z r -T- (x, -y) kM 



kM 



'^l^M^-1^^1 



(H-16) 



where 






.kM" 



1 l;>^-'S''M^[ikZ(C)+Z'(C)]e-^iFIT<^-«> 

/M"^-l 1 0 



kMd 



• /r5-?T^=lM2-l)dMds 



/(x-c) (M^-l)d^ 



. kM^ 



- [ikZ (x-d/M^-1) + Z' (x-dv^M^-1) ] e“^M‘^-l 



(H-17) 



Note that the substitution of (x) into equations . (H-15) and 



(H-16) for ij) and ij) requires one to relate x and y t>Y 

X i X X X 
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I 

M 

j 



i 



X = Y + <^ tan B 



(H-18) 



(see Figure H-1) 

These integrals were computed using the Trapexoidal Rule 
and IBM subroutine BESJ. A listing of the program follows 
in Appendix I. 
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APPENDIX I 
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